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Abstract 

The electronic transport properties of heavy-fermion systems were calcu- 
lated based on a semiphenomenological approach to the lattice non-crossing 
approximation in the limit of infinite local correlations augmented by crystal- 
field effects. Within the scope of this calculation using the linearized Boltz- 
mann theory in the relaxation time approximation the qualitative features of 
the temperature-dependent resistivity, the magnetoresistivity and the ther- 
moelectric power can be successfully reproduced; this is exemplified by a 
comparison with experimental results on CeCu2Si2 . 
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I. INTRODUCTION 



Strong electronic correlations are a crucial aspect of the physics of heavy-fermion sys- 
tems. The local Coulomb repulsion of the 4f or 5f lattice sites causes the atomic part of 
the Hamiltonian to be of non-bilinear form. Consequently the hybridization of the local 
f-electrons with the itinerant states cannot be treated in a conventional Feynman perturba- 
tion theory Several different theoretical concepts have been developed during the last 10 
to 15 years to treat this problem. These are rooted in the somewhat simpler Kondo problem 
or single-impurity problem that is now well-understood (for recent reviews see [^). If the 
f- moments are localized on a regular lattice the problem is rendered even more difficult, 
since the mechanism of coherent scattering for temperatures much below the Kondo energy 
scale Tk has to be introduced in the theory adequately in order to describe the observed 
coherence effects in the electronic transport properties. The first fully self-consistent micro- 
scopic theory developed in this context is based on the lattice version of the Non-Crossing 
Approximation (NCA) and is therefore called Lattice-NCA (LNCA) [|. The LNCA is 
able to reproduce the basic aspects of coherent scattering at low temperatures, such as a 
quadratic dependence of the resistivity on temperature. Recently the NCA was extended by 
a systematic continuation of the perturbation expansion giving rise to an improved starting 
point for the lattice approximation in the LNCA scheme [Q. 

Considering the complicated nature of the LNCA and its numerical elaboration a sim- 
plified approach including the essential aspects of the Kondo lattice seemed to be helpful. 
For that reason a semiphenomenological description of crystal field effects in Kondo lattices 
based on the LNCA was developed which is also applicable to the dynamical mean field 
theory (for reviews see [P]). This approach is able to give further insight in the respective 
infiuences of the different relevant energy scales in heavy-fermion systems on the measured 
electronic transport quantities. 

To begin with, the LNCA is briefiy reviewed starting with the formulation of the 
Anderson-lattice Hamiltonian in a generalized form based on the irreducible represent a- 
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tions of the point group of the crystaL The results are exemphfied for a Kramers ion in a 
tetragonal crystal field. The derivation of the conduction electron Green function is based 
on a semiphenomenological approximation to the local excitation spectrum. Using the re- 
lation between the conduction electron Green function and the transport relaxation time 
the electronic transport properties are calculated within the linearized Boltzmann theory. 
Finally, the calculations are complemented by a comparison with experimental results. 



The concept of the approximate treatment of the Anderson lattice for the transport 
properties was introduced by Cox and Grewe and extended to magnetotransport by 
Lorek, Anders and Grewe 0. It is based on the local approximation and uses the local T- 
matrix calculated with the NCA/LNCA describing a single scattering event of a conduction 
electron off the /-electrons. In the next section this approximate treatment will be briefly 
reviewed and completed by a consideration of the crystal field level scheme. The transport 
coefficients are calculated using the transport integral method of the linearized Boltzmann 
theory. 



The starting point of our derivation is given by the Anderson-lattice Hamiltonian. A 
spin-degenerate conduction electron band couples to localized crystal field split ionic states 
by means of hybridization matrix elements. Fluctuations on the ionic sites are limited to 
those between singly occupied or empty states, i. e. the repulsive interaction on the ionic 
sites is assumed to be infinite. The Hamiltonian is given by the following expression 



II. DESCRIPTION OF THE MODEL 



A. The LNCA including the crystal field 





(1) 
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The first part describes the kinetic energy of the uncorrelated band electrons. The energy 
of occupied local states is given by the second part. Finally, the hybridization between 
local and itinerant states is formulated by the third part of the Hamiltonian. F denotes an 
irreducible representation of the point group of the crystal, a a state of that representation 
and X^ia' Ta — |r'a"^ >< Fa^^l the Hubbard projection operator at the site v. 

In the periodic Anderson model the T-matrix Tka{z) is related to the local Green function 
{k, z) in the following way [H 



Tka{z)= VQ^ra{k(y)FT'a>,Ta{k,z)V^,^,Q{hj) 
V'a',Va 

which can be written in a more compact form by using a matrix formalism 



(2) 



Tka{z) = V!{ka)g{k,z)V{ka) 
introducing the hybridization vector V^ihj) and the Green function matrix F_{k, z) 



(3) 



V\ka) ^ {V^\^,,,{ka), Vr\^^^,{ka), . . . , Vl^^^,{ka)) 



(4) 



{N: number of lattice sites) 



With the exact equation of motion |(T2| for the conduction electron Green function 



(0), 
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Z - eka- ^ka{z) 

we obtain the following expression for the band self-energy 0: 

V\ka)Eik, z)V{ka) 



^ka{z) 



I + GflV_\hj)g{k, z)V{hj) 

OD 

Y!{ka)m, z)m<y) E(-l)' \G^ZY^{kamk, z)V{ka) 



V!{ka)m,z)Yi-K 



V{k(r)Gf^Y!{k(T) E{k, z)V{k(T) 



=A_{ka,z) 



(5) 



(6) 



(7) 



The aim of the LNCA is to obtain an approximate expression for the T-matrix and 
local Green function [^]. Unlike the NCA, in the lattice the formal delocalization of a 
local f-electron via hybridization into a band electron state may be followed by additional 
scattering events on different lattice sites before the re-localization into any local state takes 
place. In order to allow for these scattering processes the concept of the effective lattice site 
is introduced (all quantities refering to it are indicated with a tilde). The scattering is 
subdivided into local and nonlocal events, with the hybridization events starting and ending 
at a local lattice site, respectively. The summation of all pseudo-local scattering events 
results in the scattering matrix V_*^{ka)F^{z)V_{k(r) of an effective site. Accounting for the 
exclusive-site condition only for subsequent hybridization events the effective site Green 
function F_ determines the local Green function according to 

F(k, z) = ^ (8) 

with the definition 

AW = iEA(&^.^) (9) 

kcF 

Inserting equation (|]) in equation (|^) yields the following expression for the band self-energy 
that is central for the discussion of the transport properties within the scope of the linearized 
Boltzmann transport theory: 

= v\k<y) zm<y) (10) 

Exz) ^ + h.\z) — K\k — a, z) 
Note the fact that even a local approximation like the LNCA can include the leading 
anisotropies in the conduction electron self-energy by using an angular dependent hybridiza- 
tion vector Vika). 

B. Example: Three doublets; Kramers ion in tetragonal crystal field 



For further discussion we consider the level scheme of a Kramers ion (e. g. Ce^^) in 
a tetragonal crystal field. In this case the Hund ground state multiplet J = 5/2 splits 
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into three magnetic doublets. The hybridization is assumed to be spin-conserving without 
fc-dependence. Since the local pseudo-spin is a good quantum number, the hybridization 
couples only to one crystal field eigenstate in either of the three doublets for a given band 
electron spin direction. Consequently, the hybridization matrix is block diagonal and the 
6 X 6-matrix decomposes into two 3 x 3-matrices with vanishing matrix elements in the 
upper/lower block matrix for local down/up spin direction. Since the band electron Green 
function has to be diagonal in the spin the same structure is given for the matrix A(A;cr, z) 
resulting again in a block diagonal form of A,{z). Finally, the effective site Green function 
matrix is diagonal due to the local hybridization. As a consequence the band self-energy 
decomposes for a given band electron spin into a sum of the self-energy contributions for 
the three doublet states 

A crossover to the single impurity Kondo effect can be accomplished by furnishing the 
scattering matrix with a prefactor Cj^p that represents the concentration of the magnetic 
impurities 

V.2f,(^) ^ (12) 



After the expansion of the band self energy in equation (1TT|) with respect to the small 
quantity Qmp the following result is obtained 

^ka{z) = Q^p \J2 V-Fj{z)j (13) 



C. Approximate treatment of the LNCA 

The local excitation spectrum is given by the spectral functions l/vr ■ lmFj{z). The 
dominant infiuence on the transport properties is exerted by the many particle resonance 
structure near the Fermi surface. In the following this resonance is approximated by a 
Lorentzian curve 
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F,M = <nilSjL (14) 

which is in good correspondence to the results of a self-consistent solution within the PNCA 
procedure The low energy scale Tkq determines both the width 7 = tt/ {2N + 1) ■ ksTxa 
and the position rj = ksTxa of the main resonance. The exact position of the resonance 
is difficult to specify accurately since new results within the PNCA reveal a shift of the 
resonance structure towards the Fermi energy in accordance with Friedel's sum rule for a 
decreasing degeneracy Nj, unlike the former NCA result The position and width of the 
Nj degenerated excited crystal field levels are given by Aoj+T^-o and 7^ = vr/ {2Nj + l)-kBTKj 
(Aoj crystal-field splitting), respectively. 

As is well-known from renormalization theory the low temperature properties of a Kondo 
system are dominated by the low temperature energy scale Tk- The argument of the 
temperature-dependent function a{T /Txj) is therefore given by a relative temperature for 
the respective crystal field levels. All calculations performed in the following sections are 
based on the following temperature dependence 

:^{l-h{T/TK,f) : r«Tx, 



a{T/TK,) 



a (1 — In (T/Tjfj) N . y 



(15) 



The quadratic temperature dependence for T ^ T^j is motivated by the Fermi liquid 
character of the quasi-particle excitations at low temperature. The prefactor 7j/Aj (A^ = 
TTV^N{eF) Anderson width) is a consequence of the limiting behavior of the local density of 



states for T |T3| 



1 ~ 72 
Njito) = ImFAcu - i5) = — — ^-^ ^ (16) 

It refiects the many body nature of the resonance: only a fraction of < 1 of /- 

electrons participate in the scattering. In the temperature region T > T^j we use the results 
of the parquet diagram expansion of the Nagaoka-Suhl equations for the sd-model The 
parameters oq and b are fixed by adjusting equation ( [I3| ) to the temperature dependence of 
the resonance height within the PNCA procedure {Nj = 2). In the region between the low 



temperature and high temperature hmiting cases a polynomial fit based again on the PNCA 
result is used. 

The influence of a magnetic fleld on the transport properties can be accounted for by 
introducing a Kondo fleld Bkj that is given by the following relation for the respective 
crystal fleld level 

ksTRj = I^jBrj (17) 

The argument of the function a{xj) is now given by a generalized form that includes temper- 
ature and fleld on an equal base and guarantees an isotropic fleld dependence in accordance 
with the assumed fc-independent hybridization 



T V ( B V 



(18) 



Within the Lorentz approximation of the effective site Green function a Zeeman splitting of 
the crystal fleld states has to be taken into account 

Nevertheless, in the paramagnetic phase of the model we expect that there is only one 
low temperature energy scale, the lattice Kondo temperature, which we set equal to Tko- 
Equation (|1^) is only a natural way of parameterization of the decrease of the quasi-particle 
spectral weight with temperature or magnetic fleld and does not reflect a scaling law in the 



strict sense. Additional low energy scales like the Neel temperature |10] or the supercon- 



ducting Tc |Tl| can be calculated within the LNCA by analyzing the residual quasi-particle 
interactions [ffl. 



D. Relation to transport theory 

Based on the results of the preceding sections the band electron Green-function is now 
used to give explicit expressions for different electronic transport properties. Within the 
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Green-function formalism of the transport theory a comparison with the relaxation time 
approximation of the linearized Boltzmann equation for the isotropic (cubic) case gives the 
following expression for the relaxation time |I5| 



band 



[lmG'o-(a;, z)]'^dLJ 



(20) 



where the ^-dependence is averaged out. 

This integral can be solved analytically for a symmetric energy band of the width 2D 
yielding 

1 " 



h f 

T„{z = x-i6) = - Im — 

TT J-D \ X ~ to — UJ — l^cr{X — to) 



duj 



h 



uj — {x — ReSo-(x — i6)) 



+ 



27r [ [x — ReSo-(x — id) — uj]"^ + [ImSo-(x — id)]"^ 

u — {x — ReSa-(x — i6)) 



(21) 



arctan 



lmTi„{x — iS) 



D 



-D 



In the Fermi-liquid regime for T ^ the quasi-particle lifetime at the Fermi energy diverges 



resulting in lmT,fj{x — i6) and the second part in equation (21) dominates the relaxation 
time: 

n 



Tcr{z = X — id) 



(22) 



2lmS^(x — i6) 

Summing the spin-dependent relaxation times over the spin index leads to the total transport 
relaxation time 

t{z = X — i6) = ^ ro-(z = a; — i6). (23) 
This can be used to determine several transport quantities like the specific resistivity p, the 



— oo 
2^* 



1 T'^(u)(hu)''dLJ 
OUJ / 



thermopower S and the Hall coefficient Rh based on the transport integrals L„in ||T6 

Gir'^m* 1 

P = 
S = 

Rh = 
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1 Ln 
|e|TLio 
2 L20 



ne 



(24) 
(25) 
(26) 
(27) 
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To conclude this section we give an explicit expression for the matrix Ajj(^) that is 
needed to calculate the band self-energy in equation ([TT|). For a Gaussian density of states 
symmetric around the Fermi energy p(e) = plep) exp {—uj'^/D'^) the fc-sum in Ajj{z) can be 
performed as an energy integral 

A,,(z) = ^i: = -^nV}p{ep)w{z/D) (28) 



with the error function w{z/D) [|T^. Since the energy width of the Kondo resonance is given 
by the small energy scale ksTxj ^ D the approximation w(a; <^ 1) ~ 1 can safely be used; 
this leads to 

Ajj{z) ~ -mV^^piep) = -lAj (29) 

Consequently, within the three doublet scenario, the calculation of the transport quantities 
in the subsequent section is based upon the following self-energy: 

S,,...±.,).^JB^±^ (30) 



In the low temperature region equation ( |30D reflects Matthiessen's law, according to which 
the resistivity is proportional to the sum of the inverse relaxation times. 

III. COMPARISON WITH EXPERIMENTAL RESULTS 

In the following we compare the results of the model calculation with the transport 
properties of the heavy-fermion system CeCu2Si2. In CeCu2Si2 the crystal-field splitting 
between the low lying Kramer's doublet and two energetically-degenerated excited doublets 
is much larger than the low temperature energy scale Tko ||18| . Possible magnetic exchange 



scattering contributions to the resistivity are neglected. Unlike the crystal-field splittings 
Aoj and the Kondo temperature Tkq, the Kondo temperatures of the excited doublets Txj 
and the corresponding Anderson widths of the respective levels Aj are not known. In order 
to derive Anderson widths which are consistent with the fixed Kondo temperature of the 
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ground doublet and the Kondo temperatures of the excited doublets the so called "poor 
man's" scaling result [|2[ 

exp (-7r|e/|/2A) (31) 

is used with fixed values for the cutoff parameter D = 2.5 eV and the position of the localized 
f-level if = — 2eV (typical values). For simplicity we assume equal Kondo-temperatures 
Tki = Tk2 for the energetically-degenerated excited two doublets. The Kondo temperature 
for the ground doublet |0 >= —ri\ ±5/2 > — ?7^| =F 3/2 > is set to 9K according 



to reference |I9[ and r] = 0.467 ||T8[. The crystal-field splittings are known from inelastic 
neutron scattering to be Aqi = A02 = 350 K for both excited levels |1 >= \/l — ±5/2 > 
+ri\ =1= 3/2 > and |2>=|±l/2>, respectively. In cubic symmetry r] = \Jl/6 = 0.408 and 
|1 > and |2 > would form the Fg quartet. Finally, the density of states at the Fermi level 
N{eF), the absolute value of Fermi wave vector kp, the effective mass of the non-correlated 
band electrons, and the number density of charge carriers n are chosen as 1/eV-states, 
0.5/A, SrUf, {nif,: free-electron mass), and 2.5 ■ 10^^/cm^, respectively. The choice of these 
parameters is guided by typical values for d-metals so as to give a good correspondence 
with the measured resistivity. Additionally we neglected the temperature dependency of the 
Lorentzian widths 7^ in order to maintain the lowest possible set of parameters. The maxima 
due to higher crystal-field states will therefore appear more pronounced. The following 
paragraph will show that the simplifying assumption of an isotropic band structure results 
in quantitative differences between the calculated and experimentally-determined higher- 
order transport coefficients. Nevertheless, the aim of the present calculation is to give a 
quahtative correspondence. 

In figure |1] a representative temperature-dependent resistivity curve is contrasted with 
the model calculation for different Kondo temperatures of the excited doublets. In order 
to facilitate a comparison with the experimental data a residual resistivity of 20 nflcm and 
the phononic resistivity contribution of a LaCu2Si2 reference-sample |]2^ was added to the 
calculated curves. 
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FIG. 1. Comparison of the temperature dependent resistivity of CeCu2Si2 (left) with the model 
calculation (right). The inset on the right shows the calculated resistivity without residual resis- 
tivity and phonon contribution as a function of T^. The resistivity data is taken from reference 

The qualitative correspondence is best for an excited doublet Kondo temperature of 
about 120 K as far as the position of the high temperature maximum in the resistivity is 
concerned. Interesting enough, the width of the observed crystal-field transition by Gore- 



mychkin et al. []18| is about 100 K due to the strong interactions of the f-electrons with the 
conduction band in good correspondance with our estimate of the Kondo temperature of 
the excited doublets. 
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FIG. 2. Comparison of the temperature dependence of the magnetoresistivity |20|, the thermo- 
electric power[^] (CeCu2Si2) and ||2^ (Lao.9Ceo.1Cu2.2Si2), and the Hall coefficient of CeCu2Si2 
(left) with the model calculations (right). 



The low temperature maximum is more pronounced in the calculation. This is also true 

was 



for CeCu2Si2 samples with lower residual resistivity. The resistivity data of reference [20 



primarily taken since this represents a common feature of disordered Kondo lattices. With 
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increasing disorder in the system a distribution of hybridization strengths and positions of 
the local f-states is generated that cause a much broader distribution of Kondo tempera- 



tures according to the exponential dependence in equation (0). A distribution of Kondo 
temperatures might cause a significant rounding of the low temperature maximum. Due to 
electronic correlations the disorder-induced scattering rates grow faster than in uncorrelated 
systems. As a consequence, the coherent scattering part is reduced and the influence of dis- 



order cannot be accounted for by a simple increase in the residual resistivity . As shown 
in the inset of figure |l| (right) the low temperature part of the calculated resistivity follows 
a quadratic temperature dependence, as is expected in the Fermi-liquid regime. The low 



temperature resistivity of reference 0] reveals an interesting deviation from the quadratic 



behavior following a linear dependence. This might be caused by a a quantum critical point 



in CeCu2Si2 as recently discussed by Steglich et al. p8 . 

Based on the fixed parameter set with Tki = Tk2 = 120 K figure ^ shows a comparison 
of several additional transport quantities as function of temperature. The left part repre- 
sents experimental data from various references whereas the right part shows the calculated 
quantities. 

The calculated magnetoresistivity is in good correspondence with the experimental data. 
Even the shallow maximum around 30 K caused by the excited crystal-field states is re- 
produced. Quantitatively the magnetoresistivity at low temperatures reaches 800% in the 
calculation (with 20 nflcm residual resistivity) whereas the experimental value is 10 % . The 
reason for the reduced magnetoresistivity in the experimental data might also be traced 
back to the same influence of disorder on the coherent scattering part in the resistivity. 
Furthermore the calculation does not include possible band-structure effects. 

According to the calculation the Seebeck coefficient shows no sign-reversal; this is in 
contradiction to the experimental result. Furthermore the low temperature maximum does 
not appear in the experimental data. Furthermore, the absolute values in the calculation are 
sigificantly higher. This could be attributed to a disorder-induced smoothening of the sharp 
features at the Fermi level in the samples. Some general remarks seem to be appropriate here. 
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The temperature dependence of the thermopower of heavy-fermion systems shows a wide 
variety of different features. In the trivalent cerium systems a sign reversal with a negative 
low temperature minimum is a common feature in concentrated systems. Recently Kim and 
Cox discussed the possible realization of a two-channel Kondo impurity model in Ce^"*" 
predicting a large negative thermopower in a cubic crystal symmetry and non-Fermi-liquid 
signatures in thermodynamic quantities. In dilute Cerium systems, however, a negative 
thermopower has never been observed. In recent measurements on Lao.9Ceo.1Cu2. 2812 a 



clearly positive thermopower is observed which is in strikingly good qualitative agreement 
with our calculation (see single-impurity calculation of S{T) in figure |^). For temperatures 
T > Tk the lattice and the impurity thermopower calculated within our model are very sim- 
ilar since coherence no longer plays a role and the impurity concentration explicitly cancels 



in equation ^ On the other hand the thermopower sensitively measures the asymmetry 
of the scattering rate above and below the chemical potential. In our opinion the negative 
thermopower of CeCu2Si2, which is still not properly understood, points towards additional 
lattice effects not included in our model. Besides a possible influence of the band structure 
non-local quasiparticle interactions, which are not contained in the local approximation, 
cause significant renomalizations of the one-particle properties; this is part of the LNCA 
concept 0. These interactions mediate short-range antiferromagnetic fluctuations whose 
correlation length can grow with decreasing temperature. Assuming a quantum critical 
point (T/v — ^ 0) the fluctuations can account for the observed non-Fermi-liquid behavior in 



CeCu2Si2 above the superconducting [p^] . 

With increasing valence instability a crossover to an overall positive thermopower with 
a second low temperature maximum is observed (see e. g. in CeRu2Si2 p6[). One possible 
reason for the negative component of the thermopower is given by intersite spin-interactions. 
Consequently, due to the onset of real charge fluctuations on the f-sites the negative part of 
the thermopower is suppressed |^ . Spin-spin interactions are not included in the calculation 
so the qualitative feature with two maxima in the calculated thermopower might be more 
representative of systems like CeRu2Si2 . 
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We conclude this section with some remarks concerning the Hall coefficient. In this 
case the discrepancies between the calculation and the experimental data are especially pro- 
nounced. CeCu2Si2 shows a skew-scattering behavior typical of all heavy-fermion systems 
and this is not reproduced by the calculation. On the other hand, the Hall coefficient can 
vary appreciably in magnitude and sign within a set of samples of the same heavy-fermion 
system. The intricate influence of inter-site spin interactions might be responsible for the 
strong deviations of the calculations from the experimental data. Nevertheless, the calcu- 
lations predict in the Fermi-liquid regime a negative quadratic temperature dependence. 
With increasing temperature the skew-scattering part of the Hall coefficient increases and 
eventually overcompensates the negative Fermi-liquid contribution. As a result, a low- 
temperature minimum in the Hall coefficient develops. This is experimentally observed in 
several heavy-fermion systems and was discussed in more detail in |23 . 



IV. CONCLUSIONS 

The electronic transport properties of heavy-fermion systems were calculated on the ba- 
sis of a semiphenomenological approach to the lattice non-crossing approximation in the 
limit of infinite local Coulomb replusion augmented by crystal-field effects. The calculation 
is able to reproduce the qualitative features of the temperature-dependent resistivity, the 
magnetoresistivity and the thermoelectric power, as exemplified by a comparison with exper- 
imental data of CeCu2Si2 and Lao.9Ceo.1Cu2. 2812. The skew-scattering characteristic of the 
temperature dependent Hall coefficient is not reproduced. Nevertheless, the lack of agree- 
ment between the experimental data and the calculation could presumably be traced back 
to secondary effects which are not included in our model. The disorder-induced suppression 
of coherent scattering in the resistivity and inter-site spin interactions that might especially 
influence the thermoelectric power and the Hall coefficient are not taken into account. 

Within the scope of this approach an extension to include crystal-fleld excitations to 
magnetic singlet states can be easily accomplished by adding the respective resistivity or, 
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more generally, the self-energy contributions of the singlet states as described by Cor nut and 
Coqblin This permits the calculation of the transport properties of the uranium-based 
heavy-fermion systems, which tend to be in a 5f^-state [pO|] . 

A more sophisticated transport theory based on the Kubo formalism might be necces- 
sary to account for the generally anisotropic band structure and probably an anisotropic 
hybridization. This might result in a more quantitative agreement between the experimen- 
tal data and the calculated transport coefficients. However, in our opinion the consideration 
of non-local quasiparticle interactions is essential in order to account for the observed behav- 
ior of the thermopower and Hall coefficient. This is surely beyond the present approach and 
should be devoted to calculations based on a fully microscopic description of the electronic 
transport properties in heavy-fermion systems. 
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